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Ever since the inception of Transformation Optics (TO), new and exciting ideas have been proposed in the 
field of electromagnetics and the theory has been modified to work in such fields as acoustics and 
thermodynamics. The most well-known application of this theory is to cloaking, but another equally 
intriguing application of TO is the idea of an illusion device. Here, we propose a general method to 
transform electromagnetic waves between two arbitrary surfaces. This allows a flat surface to reproduce the 
scattering behaviour of a curved surface and vice versa, thereby giving rise to perfect optical illusion and 
cloaking devices, respectively. The performance of the proposed devices is simulated using thin effective 
media with engineered material properties. The scattering of the curved surface is shown to be reproduced 
by its fiat analogue (for illusions) and vice versa for cloaks. 
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To date, almost all examples of TO^"^ devices such as cloaking^"^^ and optical illusion^^"^^ were concerned 
primarily with three-dimensional free space waves and although those devices were successfully modelled, 
their practical applications have been hampered by the fundamental limits of their required material 
parameters which were lossy, dispersive and narrowband. In an attempt to reduce the complexity of the required 
materials of the devices, we will investigate TO for two-dimensional surfaces, which do not rely upon any 
symmetry. 

This Letter presents a general method to achieve illusion devices for surface waves, with the aim of utilising a 
flat, anisotropically loaded medium to perfectly recreate the total scattering characteristics of an arbitrarily, 
curved surface for all angles of incidence. For the cloaking device^^, a similar derivation to that of the illusion 
device is employed, and an anisotropically loaded, curved deformation is shown to behave as if it were a flat 
surface for all angles of incidence. The proposed illusion and cloaking devices, serve as proofs-of-concepts that 
can be adapted to other waves, e.g. acoustic. 

The underlying theory of the proposed illusion device is as follows. Starting with a flat surface, upon which 
electromagnetic waves are confined (e.g. a thin waveguide) we add an anisotropic layer, which modifies the 
effective space through which the waves propagate. If done properly, this material distribution^^"^^ will emulate 
the behaviour of a surface wave along a curved surface, thus giving rise to an illusion device. The proposed method 
is derived in the approximation of Geometrical Optics (GO). We write the electric field as an amplitude times a 
phase, E = ae'^, and assume that the phase of the wave varies much more quickly than the amplitude or the 
properties of the medium. Maxwell's equations are then recast as 



(la) 



-{VSxa)e'^ = B, 



(lb) 



VSxju~^iVSxa) + w^8'a = 0. 



(ic) 



Assuming that the medium is in an electrically thin waveguide, the fundamental mode is uniform for the cavity 
and that E is always normal to the top and bottom of the waveguide, we have e^/ = 3zz and a = az. Then VS lies 
exclusively in the x-y plane, meaning that equation (Ic) can be recast as 



det(/i) 



dS 



(2) 



where 
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Red probe line, 0^0 



Figure 2 | Probe line locations for illusion device validation. 



Ji 0 
0 1 



(3) 



as is noted in^^. In what follows it is assumed that the surface to be 
simulated is of the form z = (j(x,y)_ Now, in order to relate the 
properties of /i to the geometry of an arbitrary surface consider the 
eikonal equation in a general 2D coordinate system 



~ln 2 



(4) 



Comparing equations (2) and (4) we find the following relationship 
between the material properties in the waveguide and an equivalent 
2D metric tensor 



det{ju) 



Taking the determinant of both sides and using the fact that is a 2 : 
2 matrix one finds that 



1 



1 



det(/£) det(g) 
and the following equivalence can be established 

^ da\ ( da\ \ 



(6) 



f^xx f^xy 
f^xy f^yy 



dx J \dy 



dx) \dy 
fdoV 



(7) 



Having constructed a relationship between the constitutive para- 
meters of a flat surface, and an equivalent deformed surface, we 
proceed to test this equivalence. The first surface. Si, is simply a 
skewed Gaussian deformation (Fig. la) that has a height profile 
describe by Zi =e~^^^ and will act as the archetypal, rotationally 
asymmetric surface in this study. This surface deformation is 
approximately IXqX IOXq X SXq in the x, y and z-directions, respect- 
ively. The method of excitation employed herein was established 
^j^i2,23,24 ^j^j effectively a plane wave with the E-field polarized in 
the Z-direction with an amplitude of 1 V/m. Applying the equival- 
ence given in equation (7), we arrive at the required permittivity and 
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Figure 3 | Quantitative validation of proposed illusion device, (a) 9i = 0 
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permeability tensors (Fig. Ib-d) that will make a flat metallic surface 
coated with thin anisotropic materials appear (electrically speaking) 
to surface waves, as if it were in fact a curved isotropic surface 
described by Zi. 

In Fig. le-j, we observe the agreement between surface wave scat- 
tering (SWS) from the curved and flat surfaces (Animation 1 and 
Animation 2, respectively). As was mentioned earlier, though more 
emphasis is placed on the forward- SWS behavior with regards to 
illusion devices, the near-perfect reproduction of backward- SWS 
should not be overlooked if for no other reason, than as a perform- 
ance metric to gauge the fldelity of the proposed technique. Here it is 
also important to point out that for both illusion and cloaking devices 
the primary goal is to reproduce phase fronts that occur a number of 
wavelengths away from the surface deformation, and not necessarily 
the amplitude. In the case of the illusion device, we see that the range 
of the amplitude of is not the same for both the curved surface and 
its flattened analogue, though they are of the same order of mag- 
nitude, and in turn not inherently cause for concern. The reason for 
this difference is that in the case of the curved surface, there is 
effectively a degree of smooth-body diffraction occurring as the sur- 
face wave traverses the deformation which then leads to an interfer- 
ence pattern in the "shadow region" of the deformation. This same 
effect cannot be recreated by the flattened illusion device because it is 
borne out of GO which by definition does not take into account any 
form of diffraction. In order to truly appreciate the fidelity of the 
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Figure 4 | Multiple deformations illusion device. The method of excitation for this simulation is identical to the one described for the singular 
surface deformation (a) Isometric view of S2 (b) Syfixx- (c) Syfiyy. (d) Syfixy. Plane wave propagating (E^) from right to left {6i = 0), for curved isotropic 



medium (e) and flat anisotropic medium (f). Plane wave propagating (E^) at i 

illusion device to reproduce the forward- SWS a quantitative study 

was also conducted. Two probe lines (Fig. 2), each IS^^o in length are 

centered IS^^o away from the origin of the deformation, and oriented 

in parallel to incident surface. In the first instance, 0 is set to zero, and 

the amplitude of along the red probe line is plotted (Fig. 3a). Here 

we note excellent agreement between the curved surface and the 

illusion device. Next we change the angle of incidence of the surface 
n 

wave to - and plot the amplitude of E^ along the blue probe line 

(Fig. 2). Once again, we note an excellent level of agreement between 
the curved surface and the illusion device (Fig. 3b). Taking these 
results, and those in Fig. 1, into account one can conclude that the 



571 



for a curved isotropic medium (g) and a flat anisotropic medium (h). 



illusion device faithfully recreates the forward- SWS characteristics of 
its curved analogue. 

Next, we consider how we might apply this proposed technique to 
a real-world scenario where a perturbation is not a simple, singular 
curve, namely a surface S2, without axial or rotational symmetry. 
This surface is effectively a number of Si surfaces existing very close 
to one-another, but with no two sub-surfaces having the same height, 
origin or slope (along x or y). The net result of this juxtaposition of 
arbitrary, asymmetric, revolved Gaussian deformations is something 
akin to a group of rolling hills, and can be found in Fig. 4a. Once 
again, we note that the flat, metallic surface coated with proscribed 
material properties perfectly reproduces both the forward and back- 
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Figure 5 | Single deformation cloaking device, (a) Isometric view of Si. (b) e^jixx- (c) e,/i^. (d) Sy^^y Plane wave propagating (Ez) from right to left 

{9i = 0), for curved isotropic medium (e) and curved, loaded, anisotropic medium (f). Plane wave propagating (E^) at 9i = — for curved isotropic 
medium (g) and curved, loaded, anisotropic medium (h). 



ward-SWS characteristics of the original curved surface (Fig. 4e-h). 
It should be noted that the fact that primitive "sub-surfaces" can be 
used to build far more complex designs without sacrificing perform- 
ance bodes well for future applications. In terms of applications, an 
electromagnetic illusion device would be useful if a designer needed 
to recreate the forward- SWS characteristics of a curved surface whilst 
simultaneously trying to minimize (flatten) the space that which this 
device occupied, as might occur in the design of a surface wave 
antenna. 

Next, we derived the material parameters of a device that can 
cloak a rotationally asymmetric deformation from a surface wave. 
As one would imagine cloaking and illusion devices, for surface 
waves, rely on the same fundamental framing. We recall that for 
the illusion device we solved for the material parameters of a flat 



surface that would render it electrically deformed. Now, however, 
for the case of the cloaking device we want to solve for the mater- 
ial parameters of a deformed surface that would render it electric- 
ally flat. 

Starting from equation (Ic) and going through a derivation similar 
to that of the illusion device, we arrive at 



1 



1 + (V(7) 



/Off 



daV 
dx) 
da 
dy 



[da 



(8) 



whereas before, the surface under investigation is described as z = 
(y{xyy). 
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Red probe line, 0^0 



Figure 6 | Probe line locations for cloaking device validation. 

The test surface is once again an asymmetric Gaussian deforma- 
tion (Fig. 5a) whose profile for the sake of continuity is identical to 
that of the one used in the illusion device example. The material 
parameters to cloak such a deformation, as defined by the proposed 
solution, appear in Fig. 5b-d. As expected, all of the required relative 
permittivity and permeability values are less than unity, which allow 
the phase fronts to traverse the deformation at a faster rate than their 
fiat- space analogues, thereby giving rise to planar wave fronts in the 
forward-SWS region. The method of excitation and simulation for 
the cloak is identical to that of the illusion device, the results of which 
can be seen in Fig. 5e-h. Initially, a surface wave, is launched across 
the deformation at an angle of incidence of 0, for both the "unloaded" 
(Fig. 5e) and the "loaded" (Fig. 5f) case. True to form, the cloak 
completely renders the deformation electrically flat as is demon- 
strated by the total SWS (both forward and backward) in Fig. 5f. 
This process is repeated but for a different angle of incidence 
(Fig. 5g, h), and once again we observe a practically perfect cloaking 
performance (Animation 3). 

The minor amplitude variations of (Fig. 5h) that appear in the 
"shadow-region" are due to the underlying GO approximation 
employed in the proposed solution. Such an approximation places 
limitations on how fast the phase of a wave can change relative to the 
material parameters of the media in which the wave is propagating. 
The operating frequency is held constant through both of the afore- 
mentioned simulations, but the path at which the wave traverses now 
changes depending on the angle of incidence of the surface wave 
(owing to the rotationally asymmetric nature of the surface). 
Keeping this in mind, the minor shadowing displayed in Fig. 5h 
can be accounted for by the rate of change of the material parameters 
on the cloak exceeding the allowable limits set by GO. Keeping in line 
with the quantitative method used to validate the proposed illusion 
device's performance we do the same for the cloaking device (Fig. 6). 
Here the two probe lines, each IOXq in length, whose centers are ISXq 
away from the origin of the deformation, are oriented as is shown in 
Fig. 6. In the first instance, 0 is set to zero, and the amplitude of 
along the red probe line is plotted (Fig. 7a). Here we note that, as 
expected, due to the interference pattern created in the wake of the 
uncloaked surface deformation, the amplitude of E^ varies as it tra- 
verses down the probe-line (black curve in Fig. 7a). The cloaked 
surface deformation however, perfectly reproduces the characteris- 
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Figure 7 | Quantitative validation of proposed cloaking device, (a) 0j = 0 



tics of a plane wave (red curve in Fig. 7a), as one would expect. Next, 

n 

the angle of incidence of the surface wave is set to - and the ampli- 
tude of Ez along the blue probe line is plotted (Fig. 6). Once again, we 
see the expected interference pattern created in the wake of the 
uncloaked deformation (black curve in Fig. 7b), while the cloaked 
deformation perfectly recreates the behavior of a plane-wave (red 
curve in Fig. 7b). Taking these results, and those in Fig. 5, into 
account one can conclude that the proposed cloaking device faith- 
fully cloaks the surface deformation from surface waves. 

In this Letter we have proposed a general method to engineer 
arbitrary illusion and cloaking devices for surface waves. For the case 
of illusion devices, it is achieved by utilising a thin sub -wavelength 
coating of an anisotropic material to perfectly recreate the total scat- 
tering characteristics of an arbitrary surface for all angles of incidence 
as was demonstrated in a full -wave electromagnetic solver 
(COMSOL 4.3b). For the cloaking device, a similar derivation to that 
of the illusion device is employed, and an asymmetric deformation 
coated with a anisotropic material, is shown to behave as if were a flat 
surface for all angles of incidence. The proposed illusion and cloaking 
devices, serve as proofs-of-concept for a highly generalized technique 
that can be adapted to other waves in fields such as thermodynamics 
and acoustics. 
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